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In order to examine the influences of the spatial dependence of the pair potential, the d.c. 
Josephson current in d-wave superconductors is calculated using self-consistently determined pair 
potentials. The results show that the suppression of the d-wave pair potential near the insulator does 
not have serious effect on the properties of the Josephson current. On the other hand, drastic changes 
are obtained due to the inducement of a subdominant s-wave component, which spontaneously 
breaks time reversal symmetry. Especially, a rapid enhancement of the Josephson current at low 
£h ■ temperature predicted in previous formulas is strongly suppressed. 
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Nowadays to clarify the properties of the tunneling junction in d-wave superconductors is an intriguing problem, 
since there are several evidences which support the d-wave symmetry of the pair potentials in high-Tc superconductors 
. The observation of the 7r-j unction using Josephson junctions serves as a very convincing evidence for the d-wave 
symmetry of high-Tc superconductors 2 . Moreover, recent theories have shown the existence of a novel interference 
effect in the quasiparticle tunneling 3-5 . On the surfaces of d-wave superconductors, zero energy states (ZES) is 
formed at the surfaces as the results of sign change in the internal phases of the pair potential depending on the 
orientation of the surface. The ZES is detectable by tunneling spectroscopy as conductance peaks. By extending 
O , this concept to the general anisotropic superconductors, the phase-sensitivity of the tunneling spectroscopy has been 
predicted. The formation of the ZES at the surfaces of high-T c superconductors have been actually observed as zero 
\ energy peaks (ZEP) in several tunneling spectroscopy experiments 5_8 . 

. As for the d.c. Josephson effect in d-wave superconductors, two effects, i.e. the internal phase of the pair potential 
OS ' and the ZES formation, should be taken into account. In the previous papers, we have presented a general formula 
(referred to as TK formula) for the d.c. Josephson current, which includes both of these effects 9 for anisotropic 
superconductors. This theory is based on a microscopic basis, and the current is represented in terms of the coefficients 
of the Andreev reflection 10-12 . According to TK formula, for a fixed phase difference between two superconductors, 
qq ■ the component of the Josephson current becomes either positive or negative depending on the injection angle of 
Q\ [ the quasiparticle under the influence of the internal phase. Moreover, the maximum Josephson current Ic{T) is 
■ shown to have an anomalous temperature (T) dependence, that is, Ic(T) becomes to be proportional to T _1 at low 
temperatures when the ZES is formed on both electrodes 9,13 . These two effects peculiar to d-wave pair potential 
produce a non- monotonous temperature dependence of Ic{T) which is completely different from ordinary s-wave cases 
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^ \ to obtain more realistic results. There are mainly two effects one should take into account. One is the suppression 
• • ■ of the pair potential amplitude near the insulator. Similar effect has already been discussed based on the quasi- 
. £h ' classical Green's function method at the surface of d-wave superconductors 15 ~ 17 . The vanishing of the magnitude 
of the d-wave pair potential at the surface has been obtained when the angle between the normal to the interface 
and the crystal axis becomes 7r/4 + nn/2, n being integers. The other is a possibility of spontaneous time-reversal 
symmetry breaking near the insulator 18 . Under the existence of a weak s-wave pairing interaction, a subdominant 
s-wave component of the pair potential is induced near the interface, while the bulk symmetry remains pure d-wave 
15 ' 19 . Since the phase difference of the d-wave and s-wave components is not a multiple of ir, the mixed state breaks 
time-reversal symmetry. Under the influence of the subdominant s-wave component, the ZES splits into two levels. 
The possibility of spontaneous time-reversal symmetry breaking has been pointed out by tunneling experiments on 
the surface of YBCO 20 ' 21 . Since the anomalous temperature dependence of the d.c. Josephson current predicted 
from TK formula is responsible for the ZES formation, it is important to check the influence of the above two effects 
on the previous results. Although several properties which can not be expected from the formula based on spatially 
constant pair potentials was obtained in recent theories 22 > 23 > 24 ; the influence of the subdominant s-wave component 



The derivation of the previous TK formula was based on a spatially constant pair potential model. Although this 
model is suitable to obtain an analytical formula, the spatial dependence of the pair potential should be introduced 
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on the rapid enhancement of the Josephson current has not been clarified yet. In this paper, we will resolve this point 
in detail by extending the previous TK formula to spatially varying potential cases. 

In the following, we will calculate the Josephson current in the <i-wave superconductor / insulator / d-wave super- 
conductor (d/I/d) junction. The assumed spatial dependence of the pair potential is described as 

, „_ / A L (x,0)exp(i<p L ), (x<0) f A L (-oo,0) = A L (0) = A o cos[2(0-a)], m 
^ X ' U) ~ \ A R (x,6)cxp(i<PR), (x>0) \ A R (oo,9) = A R (9) = A cos[2(fl -/?)]. [l > 

In the above, 9 and a denotes the injection angle of the quasiparticle, and the angle between the normal to 
the interface and the crystal axis of the left (right) superconductors, respectively [see Fig. 1(b) of Ref. 9(b)]. The 
insulator which is located between the two superconductors is modeled by a 8 function. The magnitude of the 6 function 
denoted as H determines the transparency of the junction <7n- In the case of spatially constant pair potentials, with 
Al(x, 9) = Al{9), and A R (x, 9) = A R (9), the Josephson current obtained by TK formula is expressed as 25 



R N I(<p) = nR ^ k B T { y2 I m [G 12 (0-,0-,6,iu n )G 21 (0 + ,0 + ,6,iuj n )}t(e)cos6d6}. (2) 

e . . J -it/2 



In the above, t(8), [t(9) = un{9) / ' (p 2 N (9)-K 2 )] denotes the matrix element of the tunneling Hamiltonian with pn{9) = 
(m/nlcF cos9fi 2 ), where kp is the magnitude of the Fermi momentum in the d-w&ve superconductor. The quantity 
Gi2(0_, 0_, 6, iu n ) (G2i(0+, + , 6, iu n ) ) denotes the anomalous Green's function (conjugate of the anomalous Green's 
function) at the interface with the Matsubara frequency w n — 2irk B T(n + 0.5). They are given as 



k Fx h 1 + t)l,+t)l,- k Fx h 1 + r] R .+r] R .. 



sgn(w„) I A L (8±) I exp(±ia±) _ sgn(w„) | A R (6±) | exp(=F«/3±) ,., 
Vl,± = -p: , VR,± = ' ( 4 ) 

Ul n + il n ,L,± w n + "n,_R,± 



n n ,R(L),±=Bga(w n )^A 2 R[L) (6 ± )+w*, 9+ = 9, 0_ = tt - (5) 



A L {9 ± ) | A R (9 ± ) | 



The quantity Rjy denotes the normal resistance and is expressed as 

i r /2 _ OM _ . ™s 2 9 r7 _mH 



R N = / (TNCOs6d6, (XjV = — 2Q ; „ 2 , Z = — y- (7) 

Let us move to the case where the pair potential has a spatial dependence. Even if we take into account the spatial 
dependence of the pair potential, Eq. (1) does not change. Only the anomalous Green's functions Gi2(0_, 0_, 6, iu n ) 
and (G2i(0 + , + , 6, iuj n ) ) arc influenced by the spatial dependence of the pair potentials. The resulting Josephson 
current can be expressed as 

RnI ^) = *RNk B T { s^ r /2 p {ejuJnjip)aNCOsede}j (8) 

p (B ju] Lo) = 4r; L , + (0-, 0)t)r, + (O + , 9) exp(zy) 

with f = ipL — ip R . In the above, r]L,±{0-,&) and ^±(0+, 9) are obtained by solving the following equations 15 ' 16 

- ^Vl,+ (x, °) = * — ~ — a l^L ( x > + )r,l, + {x, 9) - A L (x, 9 + ) + 2w nVL , + (x, 9)} (10) 
ax Tivf cos 8 
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- £-VL,-W) = *„ l ~ [A L (x,O-)r,l_(x,0) A* L {x,6.) + 2uj n r ]L ^{x,6)\ (11) 



^VR,+ & 6) = 1 [A R (x, + )r? H>+ (x, 9) - A* R (x, 9+) + 2u; nVRt+ (x, 9)} (12) 



^Vn,R,-& 9) = hvFCOS ^ A R^ 0) - A ^ M + 2 "nVR,- 9)}. (13) 

with the boundary conditions t)l,±(— oo, 9) = fji,,± and i] R ^±(oo, 9) = f} Ri ±. In the above derivations, we have assumed 
that the spatial dependence of the pair potentials A R (x,9±) and Al(x,0±) can be determined independently as a 
semi-infinite superconductor based on the quasi-classical Green's function method 15 ' 16 ' 23 ' 24 . This assumption is valid 
for low transparent junctions (a^ <C 1). 

In the following, I((p) is calculated for symmetric junction configuration, i.e. cases where a — —ft are satisfied. The 
quantities A R (x,9±) and Al(x,0±) are decomposed as 

A R {x, 9±) = A R4 (x) cos[2(0 ± a)} + A R . s {x), A L {x, 9±) = A L4 (x) cos[2(0 =f a)] + A LjS (x). (14) 

Since the temperature dependence of the pair potential amplitude is assumed to follow ordinary BCS relation, the Tc's 
for s-wave and d-w&ve components (referred to as T s and Tj) directly correspond to the magnitude of the attractive 
interaction for the two component. We will consider two cases; (1) Al^(-i) = A* R d (x), Al iS (—x) = A* R s (x) and (2) 
Al,<i{— x) = A R ^(x), Al iS (—x) = A RiS (x) with positive x. In general, the relative phase between the s and d-wave 
components of the pair potentials, i.e., A R ^ d {x) and A R i L \ s (x), is not a multiple of tt, which implies that the time 
reversal symmetry is broken. The extreme case is a = ±7r/4, where the relative phase becomes ±7r/2. In Figs. 1 and 
2, the temperature dependence of Ic(T) is plotted for the case (1). As a reference, the results for the spatially constant 
pair potential case are plotted as curves A in each figures. The transition temperature for the s-wave pair potential T s 
without rf x 2_ y 2-wave pair potential is chosen as T s = for curves B in Figs.l and 2, T s = 0.09Td for curve C in Fig. 2, 
and T s = Q.llTd for curve C in Fig.l and curve D in Fig. 2. For these parameters, s-wave component is subdominantly 
induced only near the interface. For a = O.lw, the non- monotonous temperature dependence of Ic(T) is obtained 
even if the spatial dependence of the d-wave component of the pair potential is taken into account (see curve B in 
Fig. 1). The deviation of Ic(T) from curve A is much more drastic for a = 7r/4 (see curve B in Fig. 2) where the 
suppression of the pair potential is most significant 13 > 22 . The result shows that although the amplitude of the current 
is reduced, the qualitative feature in the Ic(T) does not change seriously as compared to curve A 13 > 22 . On the other 
hand, when the s-wave component, which breaks the time reversal symmetry, is induced, the enhancement of Ic(T) 
at low temperatures is significantly suppressed (see curve C in Figs. 1 and 2). With the increase of the magnitude 
of the induced s-wave component, the suppression of Ic{T) becomes much more significant (see curve D in Fig. 2). 
These features are intuitively explained as follows. Under the existence of the s-wave component which breaks the 
time reversal symmetry, ZEP in local density of states (LDOS) is split into two, and LDOS at zero energy is reduced. 
At the same time, the magnitude of the anomalous Green's functions Gi2(0_, 0_, 8, iu n ) and (G2i(0+, + , 6, iuj n ) ) 
for low energies decreases as compared to that for T s = 0. Thus the enhancement of Ic(T) is suppressed. 

In Fig. 3, the effect of the relative phase of the induced s-wave component between the left and the right su- 
perconductors is discussed for symmetric junction configuration (a = —/3 = 7r/4). In the case of (1), A R (x,6) and 
Al(— x, 9) are given as [Ad (x) sin (29) — iA s (x)} and [— Ad(x) sin(2#) + iA s (x)], respectively for x > 0. The relative 
phase between A R (x,9) and Al(x,6) is always ±7r independent of the magnitude of A s (x). In this case, the position 
of the free energy minima of the junction is always located at tp — ±7r independent of the temperature and Ic(T) 
increases monotonically with the decrease of the temperature (see curve B in Fig. 3). On the other hand, for case (2), 
A R (x, 9) and A^(— x, 9) are given as [Ad(x) sin(2#) + iA s (x)] and [— A d (x) sin(20) + iA s (x)], respectively for x > 0. 
The relative phase between A R (x,9) and Al(x,9) deviates from ±tt depending on the magnitude of A s (x). The 
Josephson current carried by the induced s-wave component of the pair potential flows in the direction opposite to 
that of the cZ-wave component. Consequently, the position of the free energy minima of this junction changes from tt 
to with the decrease of the temperature. This fact arouses the non- monotonous temperature dependence of Iq (T) 
(see curve C in Fig. 3). 

In this paper, we have discussed the effect of the spatial dependence of the pair potential on the d.c. Josephson 
current in d/I/d junctions. Although the suppression of the pair potential near the insulator affects the magnitude 22 
of Ic(T), the temperature dependence of Ic(T) is similar to that in spatially constant pair potential cases. On the 
other hand, the inducement of the subdominant s-wave component which breaks the time reversal symmetry has a 
large influence on the qualitative features of the Josephson effect. Especially, the enhancement of the Ic(T) at low 
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temperatures due to the ZES is strongly suppressed. Also it is a remarkable fact that the relative phase between the 
induced s-wave components of the pair potentials of the left and right superconductors can be distinguished through 
the temperature dependence of IciT). 

One of the serious problems in the field of the Josephson effect in high-T c superconductors is the lack of the 
consistency between theories and experiments. By comparing the present results with experimental data in detail, 
we expect to obtain consistency with the microscopic information of the paring interaction strength. Also it is an 
interesting problem to clarify how the roughness at the interface 13,22 influences the subdominant pair potential and 
the d.c. Josephson effect. 
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FIG. 1. Maximum Josephson current Ic(T) in d/I/d junction with a — —f} — Q.ln plotted as a function of temperature 
with Z = 10. A: Step function model, B: Self-consistently determined pair potential without s-wave component, T s = C: 
Self-consistently determined pair potential with s-wave component. T s — 0.17Xd. 

FIG. 2. Maximum Josephson current Ic(T) in d/I/d junction with a = —f3 — 0.257T plotted as a function of temperature 
with Z — 10. A: Step function model, B: Self-consistently determined pair potential without s-wave component, T s — 0, C: 
Self-consistently determined pair potential with s-wave component. T s = 0.09Td. D: Self-consistently determined pair potential 
with s-wave component. T s — 0.17Td- 

FIG. 3. Maximum Josephson current Ic{T) in d/I/d junction with a = —/3 = 0.257T plotted as a function of temperature with 
Z — 10. A: Step function model, B: Self-consistently determined pair potential under the coexistence of the s-wave component, 
with A ijd (— x) = A* R d (x), Al, s (—x) = A* R s (x). C: Self-consistently determined pair potential with s-wave component with 
A L ,d(-x) = ArXx), Al !S (-x) = A R , a (x). In B and C, T s is chosen as T s = 0.09T d . 
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